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Abstract 



Interplay of nonlinear and quantum effects in the ground state of the 
E(K)(6i + 62) Jahn- Teller model was investigated by the variational approach 
and exact numerical simulations. They result in the recognition of (i) im- 
portance of the admixture of the first excited state of the displaced harmonic 
oscillator of the symmetric phonon mode in the ground state of the system 
in the selftrapping-dominated regime; (ii) existence of the region of localized 
bi-undisplaced oscillator states in the tunneling-dominated regime. The ef- 
fect (i) occurs owing to significant decrease of the ground state energy on 
account of the overlapping contribution of the symmetric phonon mode be- 
tween the states of the same parity. This contribution considerably improves 
variational results especially in the selftrapping-dominated regime. Close to 
the E(8)e limit, the nonlinear effects of two-mode correlations turn to be effec- 
tive due to the rotational symmetry of this case. In the tunneling-dominated 
regime the phonon wave functions behave like the strongly localized harmonic 
oscillator ground state and the effect (i) looses its significance. 
pacs[63.20.Kr,31.30.Gs,71.70.Ej ] 

I. INTRODUCTION 

Recently, revival of interest in two-level electron-phonon systems occurs owing to the 
experimental evidence that Jahn- Teller structural phase transition occurs in spatially 
anisotropic complex structures (fullerides, manganite perovskites, etc)^'^'^'^'^. 

JT model is a protopype model for phonon removing the degeneracy of electron levels^' ^. 
Current investigations were focused mainly on E®e JT model with electron coupling to two 
degenerate intramolecular phonon modes, an antisymmetric and symmetric one with respect 
to the reflection. 

The reflection symmetry of two-level electron-phonon models like the exciton and the 
dimer model with onsite electron coupling to one phonon mode implies nonlinear peculiarities 
of quantum nature^. For these models. Shore et al.^ introduced variational wave function in 
a form of linear combination of the harmonic oscillator wave functions related to two levels of 
different parity with respect to the reflection. This picture can be understood in terms of two 
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or more asymmetric local minima of the effective polaron potential (i.e. the potential energy 
expression for the trial wavefunction in the space of variational parameters). Here, respective 
ground state wave function can be approximated by a linear combination of two oscillators 
with parameters corresponding to these local minima and coupled by means of further 
variational parameters. This approach was shown to yield the lowest ground state energy 
for the two-level modcls^'^. The peculiarities due to reflection phenomena occurred in the 
rotation-symmetric ECSc JT modcl^ as well. Strongly localized non-displaced phonon (exotic) 
states appeared in the numerical spectra. They were considered useful for interpretation of 
the "fast" component of luminiscence spectra. However, in order to respect the rotational 
symmetry of the E®e model the proper variational approach should be formulated in radial 
coordinates^^. 

In crystals exhibiting high spatial anisotropy with tensor properties of bulk character- 
istics (e.g. perovskites, fuUerids, etc.) the rotation symmetry of Jahn- Teller molecules is 
generally broken. Therefore, it is reasonable to investigate JT model assuming different cou- 
pling strengths a and (5 for the onsite intralevel and inter level electron-phonon couplings, 
respectively (E(8)(6i -|- 62) model^). Such a model can be also considered as a generalization 
of the exciton-phonon or the dimer-phonon model assuming the electron tunneling to be 
phonon-assisted. 

In order to understand the physical nature of the nonlinear effects we propose variational 
ansatz inspired by the shape of the numerical ground state wave function (Fig.l): the 
"principal" part of Gaussian character for both oscillators (in the absolute minimum of the 
nonlinear effective potential (Fig.l of our previous paper^^ and Fig. 2 of this paper) and a 
minor "reflective" part which corresponds to the another (local) minimum of the potential. 
For this minor part we consider the admixture of the first excited harmonic oscillator of the 
symmetric mode (rather than the only ground oscillator state, as it was commonly considered 
elsewhere). This admixture in the variation trial function leads to essential improvement of 
the results, as it will be shown in Sect. HI. 

Formulation of the variational ansatz and calculation of the ground state energy is pre- 
sented in the Section II. In the Section HI, analysis of the interplay of quantum effects 
and nonlinearity and related discussions as well as the reliability of different variational 
alternatives was investigated by comparison with results of exact numerical simulations. 



II. VARIATIONAL WAVE FUNCTION OF THE GENERALIZED JT MODEL 

We investigate local spinless double degenerate electron states linearly coupled to two 
intramolecular phonon modes described by Hamiltonian 

H = n{h\h + h\h + 1)/ + oc{h\ + h)a, - I3{h\ + h)a,, (1) 

where ar, = (^^ a-, = ( }^ ^1 are Pauh matrices, / is the unit matrix. This 

pseudospin notation refers to two-level electronic system, the Hamiltonian is thus 2x2 
matrix. 

The antisymmetric phonon mode bi splits the degenerate unperturbed electron level 
(j = 1,2) while the symmetric mode 62 mediates the electron transitions between the lev- 
els. This latter term represents phonon-assisted tunneling, a mechanism of the nonclassical 
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(nonadiabatic) nature as well as is the pure tunneling in related exciton and dimer models. 
Evidently, for (3 = 0, the one-level Holstein model (1) is harmonic, while the coupling to the 
higher level (/5 7^ 0) is the origin of the strong nonlinearity in the phonon space as will be 
seen below. 

For /3 = a, the interaction part of (1) yields the rotationally symmetric E<Sie form^ with a 
pair (an antisymmetric and a symmetric under reflection) of double degenerated vibrations. 

The general case a ^ f3 exploited here breaks the common invariance of E^e Jahn- Teller 
model under the exchange of 1 and 2 phonons; in other words, the rotational symmetry of 
the Hamiltonian (1) is broken while the reflection symmetry is kept. This symmetry is 
inherent property of the Jahn- Teller model and it is crucial for the diagonahzation of (1) 
and construction of the variational ansatz as it is evident from what follows. 

Hamiltonian (1) can be diagonalized in the electronic subspace using the Fulton- 
Gouterman unitary operator^^ 

as follows 

Hfg = UHU-' = n {bih + blb2 + 1) + a{b\ + bi)I - (5{bl + fesj^.d = H^h + Ha + H^. (3) 

The operator Gi in (2) is the phonon reflection operator: Gi{b\ + 61) = — (6| + 61), G\ = 
1. We see that Fulton-Gouterman transformation reveals high nonlinearity in the system 
(term with Gi in (3)), otherwise this nonlinearity was hidden in the initial Hamiltonian (1). 
Hamiltonian distincts from the exciton (dimer) by the phonon-2 assistance of the tunneling 
amplitude PGi. The factor Gi (2) represents continuous virtual emission and absorption 
of the phonons 1 and mediates Rabi oscillations of the electron between the levels. These 
quantum oscillations are essentially the origin of the nonlinearity of the reflection symmetric 
model as will be seen below. 

The full reflection operator is G = C^'Gi, where the electron reflection operator is defined 
by C^'ll) = |2). Equivalently to the FG transformation one can exploit commutation of G 
with Hamiltonian (1), [H, G] = 0, so that the wave function of Hamiltonian (1) related to the 
representation of the inversion group p = ±1 is a linear combination of the base functions 

|*(^)) = i=(l+pG)|l)|0(^)), (4) 

where Gi\(f)^^^) = l^*-^^-*). Hamiltonian (3) though diagonalized is no more reflection sym- 
metric, but the interaction part is antisymmetric against reflection. As a consequence, in 
the limit a = (5 the rotation symmetry is broken. Therefore, use of the FG transformation 
in the case of E(8)e JT is inappropriate: it breaks the symmetry which is necessary for the 
proper choice of the ground state in this case (see the discussion in the Conclusion). 

Inserting the representation of the wave functions (5) into the Schrodinger equation 
related to (1) we are left with the Fulton-Gouterman equation 

H'iU^^ = [Q (b\b^ + blb2 + 1) + a{b\ + 61) - pl3{b\ + b2)Gi)](t>^^^ = E^^(t>^^\ p = ±1. (5) 
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It is evident that procedure yielding the set (5) is equivalent to the FG transformation 
yielding Hamiltonian (3) in the Pauli 2x2 matrix representation. From the diagonalized 
form (5) we can see that /9-coupling breaks the degeneration of two electron levels which do 
not have the same energy. 

As it was pointed out to us, Hamiltonian (1) and, equivalently also (5) (and all the 
related quantities), are symmetric against simultaneous exchange of p ^ —p, a ^ —[3 and, 
ensuingly, in terms of parameters introduced later on, x — i^l^ ^ 1/x A* = o?' 1^ ^ 

In what follows we shall investigate variationally the ground state (p = 1) with lower 
energy only in the representation of Eq. (5). 

In two-level electron-phonon systems with linear coupling the coherent phonon subsystem 
does not conserve the number of phonons. Therefore, the upper level will share partly 
the distribution of phonons even in the ground state. From this reason, the two-center 
wave function in the form of asymmetric nonunitary ansatz with a variational parameter 
rj proposed by Shore et al^ and Sonnek et al^ for exciton or dimer models coupled to one 
phonon mode 

|*W) = -^(l + r?G)|l)|0W), (6) 

(Co is a normalization constant) was proved to yield better (lower) estimation of energy of 
the ground state''''^ when compared with the eigenfunction l^*^^^) of Hpc (5). 

Our present suggestion of the variational ansatz for the ground state is motivated by 
the numerical solution for the wave function to the diagonalized Eq. (5). In the following 
numerical and analytical calculations it is convenient to use two basic parameters: the 
asymmetry parameter x = /9/q; and the effective coupling strength /i = o? jlVt^ . The 
wavefunctions for the strong coupling n = 2 and x — 0-9) X — ^i ^^'^ X = 1-5 are given in 
Figs. la,b,c (the standard numerical simulation procedure is outlined in the next Chapter) 
where the wavefunctions are depicted in the coordinate representation in the space of two 
phonon oscillators Qi ® Q2- Three distinct forms of the phonon solutions correspond to 
small, intermediate and large values of the parameter %. 

Fig. la represents the "selftrapping" region in which contribution of the electron transi- 
tions between the levels (assisted by the oscillator 2) is small. Here the main Gaussian of the 
wave function 0"*"((5i) at negative values of Qi {p = +1) is accompanied in the Qi-subspace 
by the reflective part ri(j)^{—Qi) corresponding top = —1 parity which have already inspired 
the introduction of the nonunitary ansatz (6) representing the minor reflection respective 
to the axis Q\ = 0^. However, in Fig. la, the admixture of the first excited state of the 
oscillator 2 (coordinate Q2 ~ 0) of the parity p — —1 (displaced to the right to Qi > 0) can 
be recognized, while for p = +1 both oscillators (region Qi < 0) remain in the Gaussian 
ground state. The variational treatment which is the topic of present paper aims mostly at 
capturing the situation of Fig. la as adequate as possible. 

In Fig. lb (x = 1) which represents the case of E®e Jahn- Teller molecule, the rotational 
symmetric nature of the ground state at x = 1 is easily recognizable. We mention by passing 
that an ideal rotation symmetry of the phonon wavefunction ("mexican hat"^) which might 
be expected in the adiabatic case (that is that of big fi) of a rotationally symmetric problem 
can not be reached even for very large /x, as it was explained in detail by Eiermann et 
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al^. For the symmetric E ® e Jahn- Teller case we have merely an approximation to the 
"mexican hat" wavefunction profile spoiled by the non-zero angular momentum part which 
prevents the profile to show the complete rotational symmetry. This picture combines in it 
the features of self-trapping (Fig. la) and tunnelling (Fig. Ic) cases. The region close to 
E(8)e Jahn-Teller case appears also to show the most serious discrepancies of the variational 
treatment. 

Fig. Ic represents another limit case - that of large values of x corresponding to the 
(quantum) region of dominated phonon- assisted tunnelling. There the Gaussian form of the 
wave function is retained, although the second oscillator is displaced towards Q2 > 0, while 
phonons 1 remain undisplaced (the Gaussian is centered at Qi ~ 0). 

In Figs. 2a,b,c, we sketched three shapes of the effective potential (energy expression 
from the trial functions) controlled by the displacement parameters 71 and 72 only which 
refer to the variants of the wave functions in Figs. la,b,c, respectively. 

Thus, guided by Fig. la, the variational wave function can be proposed in the following 
nonunitary form 



4^n7i,72,n,r2, A) +7710^ ^(-7i,72,n,r2,-A) + 7720^ ^(-71, 72, n, r2, -A)l (7) 



where 7 variational parameters 7,, r^, r/j, and A are introduced and defined below. 

Phonon wave functions i = 0, 1, are supposed to be squeezed coherent and cor- 

related oscillators produced applying the set of generators on the phonon vacuum state as 
follows: 

0!,^H7i,72,ri,r2, A) = Di(±7i)^i(ri)D2(72)^2(r2)5i2(±A)|0), 
0i"H7i,72,ri,r2, A) = A(-7i)5i(ri)i^2(72)^2(r2)^i2(-A)6^|0), (8) 

(|0) is the phonon-1,2 vacuum state; indices 0, 1 at 0o,0i denote the ground and the first 
excited state of displaced phonons, respectively). 

Here we defined the generators of variational displacements 7^ 

A(7i)=exp[7,(6l-6,)], (9) 

and those of squeezings parametrized by r, 

S,{n) = exp[r,(6f - b])] (10) 

which are functions of variational parameters of displacement 7^ and squeezing for i = 1,2. 
In Eq. (5) the phonon modes 1 and 2 appear coupled in a highly non-linear way (through the 
term with /3), therefore one also includes into the ansatz (7) the mode correlation generator 

>Si2(A) = exp[A(6l6|-6i62)] (11) 

with the correlation variational parameter A. 

The functions ^o^-* and r]i(j)Q in (7-8) represent displaced and squeezed oscillators in 
Qi X Q2 space whose weight is shifted to the points (+71, +72) and (—71, +72) respectively 
(see Fig. la); The function (f)^ is merely the reflection of ^q"*"^ weighted by rji. 0o~^ = Gi4>o^^ 
(see (4-6)). 
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The function 0^"^ is the excited oscillator |0 >i |0 >2 displaced hkewise 0o~^ i^i^o ^^e 

point (—71,72) of Qi space and squeezed as well (by parameters ri,r2) with weighting 
parameter ri2. In what follows we show that introducing this admixture of the excited state 
of the oscillator-2 in (7) essentially improves variational results. 

The last variational parameter A enters in generators 512 (±A) which mix phonon modes 
together; this can be visualized as effective rotation of the two-dimensional Gaussian in the 
plane (Qi,Q2); different signs ±A keep trace of the reflection symmetry against the line 
Qi = (it is best seen from Fig. lb that left and right "hills" should rotate in the opposite 
directions) . 

The complete expression for the mean value of the Hamiltonian H — Hph + Ha + Hp (3) 
in the state (7)-(ll) is given in Appendix B. A useful representation of that expression can 
be the following decomposition which separates the contributions of "ground" and "excited" 
parts of trial functions: 

{H) = 1 [((0t + r^20^^^)i/(0o + V2<t>['^))_ 
= i [(0i//0o) + 2n2{<l>lHcl>'f^^) + r)l{<l>[-^^Hcl>[-^)] , (12) 

where 0o = 00^'' + ^i0o \ 



(13) 



The effective Hamiltonian (12) involves a highly nonlinear interplay of the variational 
parameters: the admixture of the state (f)^^ contributes by the terms due to overlapping of 
the ground and first excited state of the oscillator oc 772 and by its own excitation energy 

oc r]l 

In what follows, we investigate joint effects of quantum fluctuations and nonlinearity 
in the ground state of (12) by minimalization of the energy expression with respect to 
involved VPs 71, 72, ri, r2. A, rji, 772- Parameters of the displacement 71, 72 are defined by the 
displacement generators Di{'ji) (9), parameters of squeezing by the generators of squeezing 
Si{ri) (10), parameter of the mode correlation A by the generator of the correlation 5'i2(A) 
(11) and the parameters of asymmetry 771,772 by the linear combination (7). 



III. INTERPLAY OF QUANTUM AND NONLINEAR EFFECTS IN 
SELFTRAPPING AND TUNNELING 

As it was shown in the last Section, the reflection symmetry hidden in the original 
Hamiltonian reveals in Eq. (5) due to the diagonalization by FG transformation in a highly 
nonlinear way. This nonlinearity implies new purely quantum region of the ground state with 
strong mixing of the nonlinearity and quantum fluctuations. There are several mechanisms 
supporting the quantum (nonadiabatic) fluctuations: 

In the present model the relevant nonadiabaticity parameter is the ratio of the frequency 

2 

and the coupling parameter fl/a. The ratio of the polaron energy ^ and of the frequency 
fl, o?' 1^ = 2// is a measure of the competition between the classical (polaron selftrapping) 
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and quantum effects due to zero energy ffuctuations. Tfie quantum effects related to Q are 
tlius relevant at weak couplings /i. 

The competition between the selftrapping (a) and tunneling (P) terms results in occur- 
rence of two regions of the ground state: In the phase plane %, /i, the ground state exhibits 
two phases separated by the crossover line close to x = 1 (Figs. 3, 4, 5, see also pertaining 
discussion in our earlier paper^^). It means that the effective polaron potential exhibits 
two competing minima (Fig. 2a) governed by the model parameters ^ and x- The minima 
coincide within the border of the regions lying close to the line x = 1 (Fig. 2b). The 
phase X < 1 is selftrapping dominated, with quantum fluctuations reflected in parameters 
ri,r2,X,r]i. The phase x > 1 is the phonon 2-assisted tunneling dominated region with 
continuous virtual emission and absorption of phonons 1. This phonon- 1 exchange couples 
the levels within one minimum displaced merely by 72 due to the phonons 2. This mini- 
mum is much more sensitive to the change of model parameters /i, x as well as to quantum 
fluctuations reflected in ri,r2, A and insensitive with respect to rji, while rji ~ 0. 

From the electronic point of view, the electron in the selftrapping dominated region 
is trapped by the phonons 1, but due to the interactions mediated by the phonons 2 it 
can tunnel to the higher level. Then, owing to the reflection symmetry of the phonons 
2, continuous oscillations of the electron at simultaneous virtual emission and absorption 
of phonons 1 occur. These oscillations couple the levels and thus the electrons into pairs 
localizing them in one minimum (Fig. 2c). This mechanism was described in a recent 
paper^^ for a one- dimensional lattice model. 

An insight to the importance of different variational parameters can be gained by analyti- 
cal minimalization of Hamiltonian (12) in various approximations. Numerical considerations 
show that contributions from the quantum parameters r^, r/j, A are at least by an order smaller 
than those from the classical parameters. Including also 772, we get approximately (Q = 1): 



'1 — 7?^) 

(H) = 7? + 7i + 1 + 2^71 7-— K - 2/372£ + 2772(72^ - /3)- (14) 
Assuming all nonadiabatic parameters small and minimalizing (14) we get approximately 

71 (1 + iph^) = -a , 72 = /3£ , (15) 

where e oc exp(— 27^) (13). From these equations, approximate expressions for both regions 
(small and large x) are summarized below: 

In the "selftrapping" region a > /3 (x < 1) there is £ ^ 1 and we get: 

7i ~ -a , 72 = /3e < 1, (16) 

(that is small 72 and large negative 71); In the framework of this approximation (taking also 
r]2 small) respective ground state energy in the selftrapping region results in 



For the tunnelling dominated region, vice versa: 

a 

1 + 4/3' 



71 - -T-T^ < 1 ' £ ~ 1 - 272 , 72 = /?£ ~ /?. (18) 
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In this case the ground state energy is approximated by 



Comparing both ground state energies (17) and (19) we observe the asymmetry against the 
exchange of the a and (3 of both results due to the screening of the tunnehng term oc (3'^ by 
e — exp(— 27^) which is either vanishingly small (16) or oc 1 (18). Evidently, it is caused by 
the presence of the nonlinear Gi factor in the Pulton-Gouterman Hamiltonian (3) or (5). 

Let us examine more thoroughly the parameter 772 representing the relative weight of the 
"excited" admixture. Qualitatively important features brought by the Ansatz (7) can be 
found analytically merely from the simplified version of the Hamiltonian (12) and (37)-(40) 
with only displacements 7^ and 772. 

The equation for optimized value 772 for both regions reads exactly as 

772(1 - 4«7i) - (1 - ril)(3{l - £^) = ; (20) 
Inserting there the above expressions for 7^ (15-18) we get: 
• for the selftrapping region 

- 1 + 4a2 ' ^^^^ 



for the tunnelling region 



(77) 4:l3-ff 



71 given by (18). 

This analytical estimation shows that the admixture of the excited 2-phonons should play 
the most important part in the selftrapping region only; This conclusion is in the complete 
accordance with the shapes of the wave functions for both regions (Figs la, Ic), as well as 
with further exposed results of minimalization of variational energies (Figs. 3-5). 

Further, exploiting the influence of A and rji separately (using (12), Appendix B and (16) 
for X < 1 ), we get estimations: 

sinh(2A) ~ —2(5£^i ~ 4%// exp(— 4//), 771 ~ ^^^^ ~ x^exp(— 4/x), (23) 

where = 0^/2^1^ was defined at the beginning of the Section III. as the parameter of the 
effective interaction. (The first dependence of (23) is recognizable on Fig. 6 and will be 
discussed later in relation to the necessity of accounting for mode correlation). 

We calculated then the optimized values of the variational parameters finding numerically 
the minimum of the energy functional in different approximations - starting from the most 
complicated case of complete expression (Appendix B) with 7 parameters included and 
ending with the adiabatic ansatz with the displacements 71,72 only. 



8 



In order to check the vahdity of the variational calculations we performed also the nu- 
merical diagonalization of the Hamiltonian in the phonon-1,2 space. We truncated the 
(infinite) phonon space by Ni 1-phonon states and N2 2-phonon states, thus the state vector 
is A^i X N2 dimensional. As numerical diagonahzation results show, about 20-50 phonon 
states are sufficient for convergence. In Figs. 3 to 5 we also show the results of numerical 
diagonalization of Hamiltonian matrix as function of x for ~ 0-5, = 1 and n = 4. In the 
first two cases wc took 20x20 state vector, while in the latter case to achieve satisfactory 
convergence (especially for the tunneling- dominated region when % > 1) we had to increase 
the number of phonons-2 up to 50. 

The two energy curves from exact numerical solution and "adiabatic" variational treat- 
ment present correspondingly the lower and upper bound for variationally calculated energy, 
thus any reasonable variational results should lie between these bounds, and reliability of 
a variation ansatz for a given parameter region can be judged according to how close the 
corresponding ground state energy is to these limits (Figs. 3-5). 

We limited ourselves by showing merely some crossections of the plane x, // which seem 
to be typical in discussing the validity of the variational approach vs "true" solution via 
numerical diagonalization in the phonon subspace. From Figs. 3b to 5b showing the differ- 
ences between "variational" and "exact" energy for various variational approximations it is 
seen that while the curves with all parameters r)i,X,ri,ji included give minimal discrepancy 
(which is evident, since increasing the number of variational parameters within the same 
trial function class leads to improving the results), the maximal discrepancy, and hence the 
maximal effect of additional parameters ?7j, A is observed near the line X — this region of 
observed maximal discrepancy being shifted to the point x = 1 with growing fi. 

From the results of numerical minimalization of the variational expressions for energy (see 
Fig. 3-5) one can see that the "excitation-reflection" ansatz (ERA) (7) {r}2 included) results 
in fascinating improvement of the variational simulation of the problem in comparison with 
the "simple reflection" ansatz (SRA) (771 only); especially it is evident from Figs. 3b, 4b, 5b, 
where the differences between variational and "exact" energies arc plotted. Although this 
ansatz can be to a larger extent inspired by merely the shape of the wave function (Fig. la) 
one might wonder that the ansatz containing an excited oscillator as the "reflective" part 
of the trial function lowers the total energy of the system in comparison with the ansatz 
containing the "zero" state of the oscillator (772 = 0). An insight to better understanding 
this phenomenon can be gained examining the variational energy expression (12), (37)-(40). 
The energy can be split into three parts representing respectively the energy of the "main" 
Gaussian, that of the reflection part and "overlapping" exchange terms (those containing 
77°, r)f, r)l respectively). Indeed when we switch from SRA towards ERA the energy of the 
reflection part is increased by ~ t^I by virtue of one extra displaced "phonon" (40); But, if 
one compares the overlapping terms for both expressions (below) one can see that it is the 
overlapping term of ERA (38), (24) which significantly decreases the overall energy, while 
the respective overlapping term of SRA (37), (25) contributes only slightly. 

The main contribution to the overlapping integral is contained in the term 
(<I>o(7i)|Q2G'|<I>i(— 7i)) = ($o(7i)|Q2|*I'i(7i)) ~ 1 (Q2! second phonon coordinate, is a short- 
hand for (6^ + ^2)); other terms contain a small "overlapping" factor e ~ exp(— 27^). Using 
rough estimations (the same as those leading to (15)) we get following expression for the 
"overlapping" part of energy in the ERA: 
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El ~ -2772 (/3e^^^ cosh A + ae^'' sinh A) ~ -2r]2P , 



(24) 



(this expression is vahd for the region % < 1). In (24), it is the /5-term which yields the 
main contribution. 

As for the SRA, its counterpart reads as ($o(7i)|<52|^o(— 7i)) and vanishes due to sym- 
metry, leaving us only smaller terms oc £,: 

El ~ -Aprjie (25) 

(the higher order terms ~ e"^, e^r] etc. are omitted). 

Comparing these expressions we can see that ERA yields better use of the reflection 
symmetry property of Hamiltonian contributing greatly to the overlapping integral (to the 
negative exchange energy) while for the SRA this principal contribution vanishes merely 
because of symmetry, leaving us with minor contributions ~ e only, thus there the whole 
idea of the reflective ansatz losses much of its effectivity. 

This effect of lowering energy in the excitation ansatz due to overlapping finds its origin 
in the presence of phonon-2 assistance. In the dimer or exciton models instead of l3Q2(Jx 
term of the model Hamiltonian (1) there stands merely Actj; with the constant and 
this principal part of the exchange energy (overlapping integral) for the excitation ansatz 
vanishes. 

The very general impression from Figs 3-5 (b) makes us to state that introducing nonuni- 
tary parameters {r]i) essentially improves the variational treatment for the self-trapping re- 
gion; In the tunnelling region merely Gaussian expressions for displaced oscillators with 
squeezing (parameters rj) gives us a satisfactory fit. Indeed, as it was demonstrated from 
analytical estimations and as it is seen from Figs. 1 a,c, the admixture of reflection part is 
relevant rather for sclftrapping region where this form of the trial function is the best choice. 

However, the closer we arc to the intermediate region between sclftrapping and tunnelling 
phases, the stronger are discrepancies for all curves. As it is seen from Fig.l b it is the case 
where the wavefunctions display their radial symmetric structure. Examining on Figs. 3-5 
the curves corresponding to variational ansatzes with or without mixing parameter A we see 
that this parameter essentially lowers the energy exactly in the region of x — 1- It is worth 
comparing Fig. 3 for = 0.5 (weak coupling) with Figures 4 and 5 (yU = 1,/x = 4), both 
representing strong couplings. We see immediately that the coupling parameter A gains 
importance rather for small couplings where it improves the results for wider range of x, 
and not only for x — 1- 

Fig. 6 where we plotted the differences of the variational energy calculated with and 
without taking into account the mode correlation 5*12 (A) illustrates this statement (see also 
(20)) by showing the regions of importance of the correlation parameter A in the whole plane 
(x, /u)- The mode correlation represented by A (Fig. 6) appears to be by one order larger 
than the contribution of the competing nonlinearity due to the reflection level mixing rji. 
The correlation A is most effective for weak effective couplings /i at x — (0-5, 1.5) where it 
competes the selflocalization in support of the tunneling phase. For large 11 it contributes 
only very close to x = 1, where it reveals a maximum for all This is quite understandable 
if wc note that introducing 5*12 (A) (11) means effective rotation of trial functions (displaced 
Gaussians) in the plane Qi ® Q2, and 0+(A), (/)^(— A) arc rotated in the opposite directions 
symmetrically with respect to the line Qi = 0, which indeed repeats Fig. lb, fitting the 
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rotation symmetry features. This degree of freedom allows to represent the picture of the 
wave function especially in the transition region where selftrapping and tunnelling regions 
are mixed together and are hardly distinguishable, which is the case of weak couplings. At 
strong couplings those regions are more pronounced, the border between them is sharper; 
in this case the parameter A looses its importance with exception of the vicinity of x ~ 1 
(Fig. 6). 

In the case of omitting A, from (36)- (40) one can see that the optimized value r2 = 0, i.e. 
the contribution due to r2 in (36) is mediated merely by the correlation A. For A 7^ 0, the 
squeezing r2 significantly interplays with ri especially for small /i and x ~ 1; sts it brings 
almost half of the contribution of r2- This effect was omitted in the variational treatment 
of E(8>e model by Lo^^. Some authors disregarded this circumstance setting ri = r2 for 
simplicity, but omitting the mode correlation (A) which is therefore not selfconsistent. 

In the lattice case, the coupling with the lattice is represented by a transfer term in the 
Hamiltonian of the order of magnitude of the bandwidth T^^. When T is sufficiently large 
so that the "transfer" part of the energy is comparable to the "local" energy contribution, 
the effect of rji is considerably stronger, and we do not observe suppression caused by in- 
troducing extra correlation between phonon modes (that is an analogue to the parameter 
A). The phonon-1,2 correlation in the lattice case is of smaller order of magnitude than 
the contributions from the transfer terms (respectively of the order /?exp(— 7^) and T). Be- 
cause of that introducing the correlation VP in the lattice case does not yield considerable 
improvement of the results. 

IV. CONCLUSION 

Hamiltonian (3) allows us to distinguish two competing regimes of the electron-phonon 
system according to the relations of the parameters specifying (a) selfiocalization a vs quan- 
tum fluctuations fl and (b) tunneling P vs selfiocalization a. Then, in terms of the relevant 
parameters jj, = ^ (effective interaction) and X — ^ (asymmetry), two quantum regions 
can be identified: (i) /j < 0.5 and (ii) j3/a > 1. In these regions the quantum fluctuations 
are most pronounced and variational ansatz which pretends to be the most suitable should 
fit there the numerical data at best. Our choice of the wave function (7)-(ll) covers both 
regions in a complementary way: while the choice of ERA with the admixture of the excited 
state of the symmetric phonon mode weighted by 772 and including the mode correlation A 
improves greatly the variational results in the selftrapping region, % < 1, in the whole range 
of fl (Figs. 2, 4, and 5), with increasing fi effectiveness of all quantum variational param- 
eters vanishes except of 772 (Figs. 4 and 5). In the tunneling region, x > 1, the choice of 
ERA looses its justification in benefit of SRA, while all the remaining parameters keep their 
effectiveness even for large // (Figs. 4 and 5). 

The cooperative effect of the reflection (antisymmetric phonon mode) and of the as- 
sistance of the symmetric mode in the tunneling results in a nonlinear interplay of both 
modes. It consists in the competition between the negative contribution of the overlapping 
of the wave functions of different parity with respect to the reflection and the increase of 
excitation energy of the respective mode. This concept leads to the effective energy decrease 
of the excited symmetric reflected mode (ERA) rather than of its ground state (SRA). 
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The complex nonlinear interplay of the modes was ehicidated by exact numerical diago- 
nalization of Hamiltonian (3): in fact, we took inspiration for ERA from the exact ground 
state wave function for a > /3 (Fig. la). It suggested presence of admixture of the first 
excited oscillator state of the symmetric mode in the reflected part of the wave function. In 
the case of a < /3, the numerical wave function exhibited only a single symmetric peak of a 
well defined harmonic oscillator. The peak was located close to the center of the reflection 
symmetry Qi = but displaced by the phonons 2. It corresponded to a new minimum of 
the effective potential which opened due to the "bond selfiocalization" on account of the joint 
effect of both modes (Fig. 2c). Note that these states are well locahzed, in the contrast to 
the states in the "selftrapping" region. It is interesting to mention in this context a special 
class of states in the excited spectra of J-T models which were called "exotic states" ^'^ and 
were characterized by a pronounced localization of the corresponding phonon mode. These 
states in the excitation spectra can be explained as the consequence of the energy resonance 
and hence the tunnelling between two wells of the effective potential (visualized, e.g., on Fig. 
2) which suggested opening of the additional potential well in the position 71 ~ where the 
exotic mode is to be localized. Although in the present paper we investigated merely the 
ground state of the model, we can see that our localized modes in the tunnelling dominated 
region (x> ^) have essentially the same origin (the "localized" minimum at zero displace- 
ment along Qi of the effective potential) and bear interesting resemblance to the Wagner's 
exotic states. In support of this statement the full spectrum of the phonon states, and not 
the ground one only should be examined, but this problem merits a special paper. We just 
mention now our own results on numerical calculation of the energy spectrum for the excited 
modes in the tunnelling region which show specific periodic (in model parameter /i) chaotic 
"windows" especially for higher modes. This periodicity in the coupling strength, which 
we at the very beginning scaled by the phonon frequency clearly indicated the resonance 
behaviour, i.e. occasional coincidence of two incommensurable characteristic frequencies of 
the system whose nature can be identified with the origin of Wagner's resonant exotic states. 
For low lying modes like those of the ground state this behaviour is not so pronounced, but 
merely its traces are also recognizable. 

The results exposed above, especially when speaking about the validity of the variational 
approach chosen, are most adequate far from the rotationally symmetric Jahn- Teller case 
(a = or X = 1) which was to be expected on basis of the comment to the FG transforma- 
tion in the Section II. In this case one gets an almost degenerated degree of freedom in the 
space of variational parameters; namely, if we introduce an analogue of the polar coordinates 
in the phonon-1,2 space (as some authors^ ■^'^'^ do), certain degeneracy of the energy pro- 
file over the angular coordinate would be observed and thus this angular parameter should 
have been excluded from being variated. (Strictly speaking, the angularly degenerated true 
"mexican hat" appears only in adiabatic approximation"'). The authors exploiting various 
variational treatments disregarding this circumstance must have encountered such problems 
for the E(8)e Jahn-Teller case. However, these inconsistencies become crucial rather at big 
coupling strength e.g.. Fig. 4, 5). However, the rotational symmetry can be to some 
extent retained even within the formalism of rectangular Qi ® Q2 space by introducing 
the mode mixing parameter A which gains its significance especially in the region of weak 
coupling effectively spanning both regions (Fig. 6). The variational approach exploiting es- 
sentially the rotational symmetry of the model should present a complementary description 
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suitable in the vicinity of x = 1, the whole problem however being a subject of our further 

considerations. 
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APPENDIX A. 



We have used following formulas for D^, Si, and 5'i2 defined by (9)- (11) which can be 
found elsewhere^^'"^^: 



Di{ji)-\Di{-fi) = k + ji, i = 1, 2, 
-Di(7i)S'i(ri) = Si{ri)Di{^i), 71 = 7ie"^'', 
S^\ri)biSi{ri) = hi cosh2ri + h\ sinh2ri, 
S'{2^(A)6i5'i2(A) = hi cosh A — h^ sinh A, 

(0|5l(ri)i^I(7i)exp(A6l)) 



(cosh2ri)V2 



X 



X exp 



7? 



— tanh(2r) - A7i(tanh2r - 1) + 4j-(tanh2r - 1) 



(0|5i(r)Di(7i)6r) 



■(0|5i(ri)A(7i)exp(A6l))| 



A=0; 



T(<5) = exp [S{hlh2 - 6i4) 



Tt 



J_fhi-h2 

V2\hi + h2. 

rt|o) = r|o) = |o). 



(26) 
(27) 
(28) 
(29) 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



APPENDIX B. 



{H) 



C 



(-)t; 



;(-)\ 



(36) 



where 



{(^^o"^^ + Vi^ir^^W^^o"^ + Vi<-^) = (1 + Vl + 2r7i£) 



; (cosh 4ri + cosh 4r2) cosh 2A + 7^ + 72 



+2riie {— tanh 2A sinh 2A cosh 2(ri + r2) cosh 2(ri — r2) 
+^2^(e2{n+r2) _ g(-2(ri+r2)) cosh4A)(l + tanh=^ 2A) + 2e(-2(^i+^^)' sinh 4A tanh 2A] cosh2(ri - r2) 



cosh^ 2A 



(sinh4ri - sinh4r2) - 27i(tanh(2A)e^'^''2 '"^^72 + 71) 



+ ^(1 - riDli - ^ [(1 + ^i)[72 - 71 tanh2Ae2('-^-'-i)]£ + 27^172 
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(^J+^/f^i-)) = (cosh 4ri + cosh 4r2) cosh 2A + 7i' + 72'! + (^2...)) 

+ ^-= ((cosh 2(A + 2ri) + cosh 2(A + 2r2))] (62-) + (cosh2(A - 2ri) +cosh2(A - 2r2))(6i...)) 

H — ^ [(-sinh2(A - 2ri) - sinh2(A - 2r2))(6?...) + (sinh2(A + 2ri) + sinh 2(A + 2r2))(62--) 
4v2 

1 1 
+[-sinh4ri + sinh4r2 - -(sinh2(A - 2ri) +sinh2(A - 2r2))] — ^{b^...) 

2 2v 2 

+ [- smh4ri + sinh4r2 + ^(sinh2(A + 2ri) - sinh2(A + 2r2))]^^(6i62---) 

2 2v2 

+ [{^16^'' sinh A + 72e2'-2 cosh A)] ((6162...) + e) 
+l{-lie'^'+l2e''')e-\bl....) + ^{j,e'^^ + j,e'^^)e\bl....) - ^coshAe^'"^ 

r 2ri 

+ ^ -e^-i sinh A(£ + (6162...)) + ^ - {bl...)e') - V2^^ {{b,...) + {b^...)) 

{¥q-^^h¥-^) = (sinh A7ie2'-i + cosh A72e2'-2) - ^sinhAe^'-i - ^ [V272(< 62- > + < h... >) 
+ l(e-^+2^2 < _ ^ +e^+2''2 < 52... > +coshAe2'-2 < 6162- > +coshAe2"2 < ... > 

(vl>["^"^ifvl>(")) = ^2 _^ ^2 _^ cosh2A(cosh4ri + cosh4r2) + sinh2(ri + rg) sinh2(r2 - 

*^-27i - ^ |272Cosh2A£ + 72(< bj... > - < bj... >) sinh A 



(38) 



(39) 



,2r2 



+ 



[< bi... > e~^+ < b2... > e^)4coshA 



2V2 

+((< bl... >-< blb2... >)e^ + (< blbi... >-< bl... >)e-^)sinhA]} (40) 
The mean values in (13) 

{b^D,iV2^,e')S,{X)){D2i-V2^,e-')S2{-X)) = ^^^e; (41) 
{blD,{^/2^,e^)S,{X)){D2{-^/2^le-^)S2{-X)) 

= — ^— ftanh 2A + 27i^e^^(tanh 2A - 1)^) e; (42) 
cosh2A ^ ^ 

{blD^{V2^,e^)S,{X)){D2{-V2^ie-^)S2{-X)) 

2y,<^^it^aA^2X-ll ^^^^^^^^ ^ 272e2^(tanh2A - 1)^) e; (43) 
cosh2A ^ 

= /. = 1,2,3. 
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FIG. 2. Shapes of the effective potential corresponding to the wave functions at Figs.la, lb 
and Ic, respectively. 
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(r7i=A=0) 












(A=0) 









FIG. 3. (a) The ground state energies (36) for n = 0.5. The seltrapping dominated GS spans 
over X < 1 f-iid the tunneling dominated GS over x > 1- The curves plotted represent cases (from 
below): numerical simulation GS, Eex; ERA, Ei; ERA, Ei{X = 0); SRA, Eq; SRA, £;o(r/i = A = 0); 
adiabatic GS Ead- 

(b) Differences of the ground states from (a) and the exact numerical GS. ERA considerably 
improves the results for % < 1. It also shifts the maximum of the differences to the point X = 1. 
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b) iu=l 










(f)l=A=0 








(A=0) 







FIG. 4. The same as in Fig. 3 for /x = 1. With increasing /x, for x < 1, the loss of efficiency of 
ah VP except of r]2 is evident. 
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. 

FIG. 6. The difference of two GS: the GS without the reflection and correlation effects 
Eo{r] = 0, X = 0) and the GS Eq including these effects. 
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Figure Captions 

Fig. 1. The numerical ground state wave functions ai fi — 2 and x — 0-9 (a) x = 1 (b) 
and X = 1-5 (c). 

Fig. 2. Shapes of the effective potential corresponding to the wave functions at Figs, la, 
lb and Ic, respectively. 

Fig. 3. (a) The ground state energies (36) for fj, — 0.5. The seltrapping dominated GS 
spans over x < 1 ^ind the tunneling dominated GS over x > 1- The curves plotted represent 
cases (from below): numerical simulation GS, E^x', ERA, Ei, ERA, Ei{X = 0); SRA, Eq; 
SRA, Eo{rii = A = 0); adiabatic GS E^d- 

(b) Differences of the ground states from (a) and the exact numerical GS. ERA considerably 
improves the results for x < 1- It also shifts the maximum of the differences to the point 
X = l- 

Fig. 4. The same as in Fig. 3 for = 1. With increasing /i, for x < 1, the loss of 
efficiency of all VP except of ri2 is evident. 

Fig. 5. The same as in Fig. 3 and 4 for // = 4. 

Fig. 6. The difference of two GS: the GS without the reflection and correlation effects 
-E'o(^ = 0, A = 0) and the GS Eq including these effects. 
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